Fluctuation-dissipation ratios in the dynamics of self-assembly 
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We consider two seemingly very different self-assembly processes: formation of viral capsids, and 
crystallization of sticky discs. At low temperatures, assembly is ineffective, since there are many 
metastable disordered states, which are a source of kinetic frustration. We use fluctuation-dissipation 
ratios to extract information about the degree of this frustration. We show that our analysis is a 
useful indicator of the long term fate of the system, based on the early stages of assembly. 
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I. INTRODUCTION 

Self-assembly processes can be loosely defined as those 
in which simple building blocks assemble spontaneously 
into highly ordered structures. Assembly is of vital im- 
portance in biology, where cells use dynamically assem- 
bled protein structures to control the shapes of lipid 
membranes 

H H H, and the life cyclc of 

viruses involves 

spontaneous assembly of their protein coats 0, [f| . Re- 
cently, self-assembly has also been used to develop nanos- 
tructured materials [|| 0] , which often draw inspiration 
from biological systems. The self-assembly of viral cap- 
sids [|| has been the subject of elegant experimental Q 
and theoretical studies [a ISl lid] • In this article, we study 
a model [9( designed to mimic this assembly process. At 
low temperatures, assembly is frustrated by the presence 
of long-lived disordered states. The avoidance of this 
frustration is crucial for successful assembly. This ef- 
fect is rather general, as we illustrate by also considering 
the formation of ordered structures in a two-dimensional 
system of sticky discs. We analyze the crossover be- 
tween frustrated and unfrustrated regimes, and show that 
fluctuation-dissipation ratios (FDRs) [ll], EH E!> E3] as- 
sociated with the early stages of assembly are correlated 
with the long-time yield of these processes. This repre- 
sents a new application of FDRs, which have been dis- 
cussed previously in the context of glassy dynamics. We 
discuss how this analysis might be helpful in the design 
of self-assembling systems. 

In general, successful self-assembly requires both that 
the equilibrium state of the system is an ordered struc- 
ture, and that the system reaches this ordered state in 
the time available for the biological or experimental ap- 
plication. The first condition is thermodynamic: for the 
systems studied here, the low energy ordered states are 
known, and this crossover can be estimated by free energy 
arguments, as in @. (Note however, that if 'liquid- like' 
states are relevant near the thermodynamic crossover, 
then this will lead to more complicated behavior, as in 
fiol].) The second condition for successful assembly is 
kinetic in origin: it is illustrated for a model system 
of viral capsid assembly in Fig. [TJ The degree of as- 
sembly shows a maximum at a finite temperature T*. 
As the temperature is lowered through T*, the ordered 
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FIG. 1: (Color online) Assembly of model capsids in the B4 
model of Ref. [jj. (a) A well- formed model capsid, with icosa- 
hedral symmetry, (b) Representative selection of metastable 
states formed at reduced temperature T = 0.067 and reduced 
time t = 3x 10 s (see the text for definitions), (c) Plot of the 
the capsid yield at t = 3 x I0 5 , which is non-monotonic in the 
reduced temperature. The yield is the fraction of particles in 
complete capsids, identified as in Here and throughout, 
red and blue symbols identify high and low temperatures re- 
spectively. We also indicate the approximate location of the 
kinetic crossover, at reduced temperature T* . 



state becomes more probable at equilibrium, but the self- 
assembly process becomes less and less effective: we refer 
to this change as a "kinetic crossover" . 

The purpose of this article is to use dynamical ob- 
servables to study the behavior near T*. Since this 
is the regime of most efficient assembly, it is relevant 
both biologically and for applications of self-assembly 
in nanoscience. While the kinetic crossover can always 
be identified by measuring the degree of assembly, as in 
Fig. [TJ achieving this in a computer simulation requires 
access to long timescales, which restricts the range of sys- 
tems that can be studied. In this article, we show how 
FDRs can be used to locate the kinetic crossover using 
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FIG. 2: (Color online) Assembly of sticky discs, (a) Typical part-assembled structures at reduced time t = 5 x 10 6 , and reduced 
temperatures T = 0.17,0.2,0.26,0.27,0.33, from left to right. Illustrated regions are of size 25ao x 25oo- The crystallinity is 
poor at low temperatures, due to the metastability of the disordered states, (b) Plot of the fraction of particles with 6 bonds, 
which is a measure of the yield of the assembly process. Compare Fig. QJ. 



simulations on relatively short time scales. (It is nec- 
essary to average over many such short simulations, but 
such averaging is trivially parallelizable.) We also discuss 
how these response functions might be measured exper- 
imentally in ordering processes that occur on complex 
energy landscapes. 



II. MODELS 

A. Model capsids 

The first model that we discuss describes the assem- 
bly of viral capsids. Full details are given in The 
model consists of rigid subunits, the "capsomers" , which 
interact by isotropic repulsive forces, and directional at- 
tractions. The low energy states in the model contain 
"capsids" , each of which consists of 60 subunits in a cage 
structure, with icosahedral symmetry. We use the B4 
variant of this model, which means that the attractive 
potential favors the capsid structure shown in Fig. [TJi. 
The subunit diameter is a, and the density of subunits is 
p. The parameters of the model are the reduced capsomer 
density pa 3 and the reduced temperature T (measured in 
units of Eb/ks, where £b is the energy associated with the 
attractive potential and fcs is Boltzmann's constant). In 
addition, the specificity of the directional attractions is 
controlled by the angular parameters 9 m and <p m . The 
data of this article are obtained under the representa- 
tive conditions pa 3 = 0.11, 6 m = 1.5 and <p m — 3.14. 
We simulate a system of 1000 capsomers in a cubic box 
with periodic boundaries. The capsomers evolve accord- 
ing to overdamped Brownian dynamics, and the unit of 
time is (a 2 /48D), where D is the capsomer diffusion con- 
stant. The rotational diffusion constant of each capsomer 
is D r = 2.5(D/a 2 ), as in [|. 



B. Sticky discs 

We also consider a second model whose subunits are 
sticky discs which interact by an attractive square-well 
potential of depth J and range ao, and a repulsive hard 



core of range 0.9ao- We quench the system into the solid- 
vapor phase coexistence regime, so that the equilibrium 
state has most of the discs in a single close-packed crys- 
tallite. However, we use Monte Carlo dynamics that are 
chosen to accentuate the effects of kinetic frustration. 
We move bonded clusters as rigid bodies, allowing trans- 
lation and rotation, but no internal rearrangements. To 
reflect the slow motion of large clusters, we use an aver- 
age translational step size of 0.1(ao/M) and a rotational 
step of 7r/(10.f) radians, where M is the number of par- 
ticles in the cluster and Jap its moment of inertia (in 
units of the disc mass). The reduced time t is measured 
in Monte Carlo sweeps, and the reduced temperature T 
is measured in units of J/fcs- Clusters can rearrange 
only by bond breaking. These events are sampled by 
the 'cleaving algorithm' of [Hj], with equal fictitious and 
physical temperatures. It is an off-lattice generalization 
of the Swendsen-Wang algorithm [lfij]. At low temper- 
atures, the dynamics mean that bonds are broken very 
rarely, and aggregation of the discs is diffusion limited. 
At T = 0, the system reduces to diffusion-limited cluster 
aggregation (DLCA) [l7| . 

The crossover from effective to ineffective assembly in 
the capsid system was shown in Fig. Q] We show similar- 
results for the disc system in Fig. [51 The system contains 
400 discs in a periodic square box of side 100ao. The 
system does not reach full phase separation into close- 
packed structures on the time scales accessible to our 
simulation, so all of our data is in the out-of-equilibrium 
regime. 



III. FLUCTUATION-DISSIPATION RATIOS 

The non-monotonic yields shown in Figs. [TJ: and [2h 
mean that for the observation times considered, and 
when the temperature is small, reducing the temperature 
does not result in a decrease in of the total energy. This 
kind of 'negative response' to temperature perturbations 
is familiar in systems with activated dynamics [l4T |. In 
the self-assembling systems considered here, the kinetic 
crossover at T* is associated with a change from positive 
to negative response on the long time scales considered 
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in Figs. [TJ; and [2}d. In this section, we use fluctuation- 
dissipation ratios (FDRs) to show that the crossover be- 
tween positive and negative response has signatures that 
can be detected on much shorter time scales. 



A. Basic idea 

Fluctuation-dissipation ratios (sometimes also called 
correlation-response ratios) have been widely studied in 
the context of aging of glassy systems [ll|. Imagine ap- 
plying an instantaneous perturbation to a single subunit 
(disc or capsomer) at a time t w , and measuring the effect 
of this perturbation at some later time t. For a sys- 
tem at equilibrium, the fluctuation-dissipation theorem 
(FDT) relates the response to small perturbations to the 
relaxation of spontaneous fluctuations [H[ • 

In general, we can measure the fluctuations and re- 
sponses of any observable. Here, we focus on the the 
energy of a given subunit. In both of our model systems, 
the total energy comes from interactions between pairs 
of subunits, E = (1/2) J^. Eij, where the primed sum 
excludes terms with i = j. We denote the energy of the 
ith monomer by 



(1/2) Y, E ; 



(1) 



We measure the responses in the system as follows. 
Starting from a given initial state, the system assembles 
for a 'waiting time' t w . We then turn on a perturbation 
to the energy: SE(t) = J2i hiEi&(t — t w ), where hi is the 
(small) field applied to the ith subunit, and <d(x) is the 
unit-step function. We measure the integrated response 
to this field, 



x(t,t v 



V d((3hi) 



(2) 



h=0 



where the notation h — (hi, h%, . . .), and fi~ x is the tem- 
perature multiplied by Boltzmann's constant. The av- 
erage is over trajectories of the system in the presence 
of the perturbation. In practice, we evaluate the partial 
derivative by assigning hi — Sh to half of the subunits 
(selected at random), and hi = —5h to the other half. In 
the linear response regime (small Sh), the mean energy 
at h = can then be estimated by E(t) — N^ 1 J2i Ei(t), 
and the response by ^J-Ei(i) — E(t)]/hi. These quan- 
tities are then averaged over many independent runs of 
the dynamics. Our results for the capsid system were 
obtained at Sh = 0.05 and those for the disc system were 
obtained at (Sh/T) — 0.3. These values are small enough 
that our estimates of x(t, t w ) change very little if Sh is re- 
duced, which indicates that we are in the linear response 
regime. For systems with Monte Carlo dynamics, such 
as the disc system, the derivative in Eq. ([2]) can also be 
evaluated as a correlation function for the unperturbed 
dynamics, in which case it is no longer necessary to apply 
the field hi directly [HHJ. 
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FIG. 3: (Color online) Sample trajectories in the capsid sys- 
tem, showing ri6o (t), defined as the fraction of particles in 
bonded clusters of size 60. We use a logarithmic scale for the 
reduced time t. The fraction neo(t) reflects the number of 
capsids in the system, since disordered clusters containing ex- 
actly 60 subunits are rare. The first capsids appear at times 
around 10 4 . The system is away from global equilibrium until 
reduced times are at least of the order of 10 J . The arrow in- 
dicates the maximal time associated with our measurements 
of correlation and response functions (Figs. [4][6]) . 



For a s yste m at equilibrium, the fluctuation-dissipation 
theorem [18[ states that 



X (t,U) = C{t,t) - C(t,t w ) 
for all t and t w , where 



(3) 



C(t,t w ) = (Ei(t)Ei(t w )) - (Ei(t)){Ei(t w )), (4) 

Alternatively, we can define the response to an instanta- 
neous perturbation (impulse response), as a derivative of 
the integrated response: R(t,t w ) = —dx{t,t w )/dt w . In 
that case, the FDT states that 



R(t, t w ) 



dC(t,t w ) 
dt w 



(5) 



Away from equilibrium, we define the correlation- 
response ratio X(t,t w ) by 



R(t, t w ) — X(t, t w ) 



dC(t,t w ) 
dt w 



(6) 



Thus, X(t, t w ) is the response of the system to an instan- 
taneous perturbation, normalized by the response of an 
equilibrium system with the same fluctuations. 

In equilibrium, the fluctuation-dissipation theorem im- 
plies that X(t,t w ) = 1 for all t and t w . Away from equi- 
librium, X{t,t w ) may take any value. It is most conve- 
niently obtained from the gradient of a parametric plot of 
x(t,tw) against C(t,t w ), where the parametric variable 
is the waiting time t w [21(. We will see that paramet- 
ric plots distinguish between systems above the kinetic 
crossover region, and those below it. This application of 
the FDR is the main result of this article. 
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FIG. 4: (Color online) (a) A time-line illustrating the simu- 
lation protocol used to measure the response, (b) Response 
in the capsid system (in units of e£) at reduced temperature 
T = 0.091, as a function of time t, for i w = 0, 480, 960, 1440. 
The data are plotted with lines, since each simulation yields 
data points for all t. (c) Plot of the response as a function of 
waiting time t w , for t = 960, 1440, 1920. This is a replot of 
some of the data of the middle panel, but it allows estimation 
of the impulse response dx(t,tw)/dt w . In this case, the data 
are shown as points (squares), and points with the same value 
of i are connected by lines. 



B. Results 

In Fig. [31 we illustrate the time scale associated with 
capsid formation. The first capsids form in the system at 
times around 10 4 , and all systems shown are significantly 
out of equilibrium until reduced times at least as large 
as 10 5 . The yield measurements of Fig. [1] were taken at 
t = 3 x 10 5 . As time proceeds, the system evolves in- 
creasingly slowly towards the equilibrium state. We will 
show correlation and response data at times of order 10 3 , 
so the system is still well away from global equilibrium in 
all cases. However, we will find that systems at tempera- 
tures above the kinetic crossover region have responses in 
accordance with FDT, while those below it do not. In the 
disc system, the snapshots of Fig. show that the system 
is well away from equilibrium at times around 5 x 10 6 . 
For that system, we will show correlation-response data 
at much earlier times, those less than 10 5 . 

Some results for the capsid response function are 
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FIG. 5: (Color online) (a) A time-line indicating the simula- 
tion protocol used to measure the correlation, (b) Correlation 
function in the capsid system (in units of e£) at T = 0.091, as 
a function of reduced time t, for t w = 480. (c) Correlation as 
a function of i w , for t = 1920. The absence of time-reversal 
symmetry is clear. 

shown in Fig.[4j where we show how multiple simulations 
are used to plot the response as a function of t w for fixed 
t, which is useful for estimating the impulse response. A 
typical correlation function is shown in Fig. [5J 

Results for the FDR in both capsid and disc systems 
are shown in Fig. [SI where we have normalized both corre- 
lation and response by the equal time fluctuation C(t, t). 
[The function C(t,t) is independent of t w , so the gradi- 
ent of the parametric plot is —X(t, i w ), and is unaffected 
by the normalization.] Above the kinetic crossover, as- 
sembly is taking place, but the energy response is in ac- 
cordance with FDT, so X(t,t w ) k 1, at least for the 
times that we considered. As the system passes through 
the kinetic crossover, the FDR shrinks. While it can 
be convenient to characterize this crossover by the single 
temperature T* , it is more accurate to think of a temper- 
ature range over which the long-time behaviour of system 
changes smoothly from effective to ineffective assembly. 
This smooth change is accompanied by a smooth change 
in the FDR. 

We conclude that if a system is to be designed so that 
it assembles effectively, the correlation-response ratio can 
be used to obtain a general prediction for the regime 
of good assembly, before running the long simulations 
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FIG. 6: (Color online) Correlation response plots for (a) the 
capsid system at t — 1920, and 960 < t w < £, and (b) the 
disc system at t = 8 x 10 4 . These systems are all well away 
from equilibrium, but the response is in accordance with the 
prediction of FDT at the higher temperatures. The response 
decreases rapidly as the system passes through the kinetic 
crossover and falls out of equilibrium. The red and blue col- 
oring is consistent with that of Figs. [T]and[5] 



required to test the yield directly. 

Finally, note that we constructed Fig. [5] using data at 
constant t and variable i w , since the gradient of this plot 
gives the FDR. This procedure requires a separate sim- 
ulation for each data point. However, if we only wish 
to know if the integrated response is small or large com- 
pared to the FDT prediction, it is sufficient to use data at 
a single t w : a simple comparison of C(t, t) — C(t, t w ) and 
x{t,t w ) is already quite informative in that case (note 
however [HI). 



IV. DISCUSSION 

We now consider the kinetic and thermodynamic 
crossovers in a little more detail. We then discuss how 
the change in FDR at the kinetic crossover arises, and 
the extent to which we expect it to generalize to other 
self- assembling systems. 



A. Thermodynamic and kinetic crossovers 

We measure the yield of our assembly processes by 
running long simulations of length t y i e id (recall Figs. 1 
and 2). These simulations have three types of final state. 
At high temperatures, no assembly takes place, and the 
system consists primarily of free subunits. At low tem- 
peratures, the system evolves into a state that consists 
primarily of disordered metastable clusters. 

We also find an intermediate temperature regime, in 
which the final state has a substantial quantity of as- 
sembled products. This regime is delineated by two 
crossovers. For an operational definition of the high tem- 
perature crossover, we impose a threshold on the relative 
probabilities of bonded and free subunits at time t y i C id- 
While this definition depends on t y i e id> the position of the 
crossover has a well-defined limit as £ y i e id —> oo, which 
can be evaluated from the contribution of free subunits 
to the thermodynamic partition function of the system. 
Thus we refer to this crossover as "thermodynamic" . 

To define the low temperature crossover, we con- 
sider the relative probabilities of disordered clusters and 
correctly- assembled products at iyidd- As the tempera- 
ture is reduced, the maximum of the yield occurs when 
the disordered clusters are common enough to signifi- 
cantly impede assembly. We therefore identify this maxi- 
mum with the low temperature crossover. If we calculate 
the yield in the equilibrium state, we expect it to depend 
monotonically on the temperature, since the correctly- 
assembled states minimise the total energy in both of our 
systems. Thus, the presence of the maximum in the yield 
is a kinetic effect, that arises from the slow annealing of 
disordered clusters. This motivates our use of the term 
"kinetic crossover" . Clearly, the existence of a regime of 
efficient assembly requires that the kinetic crossover is 
not too close to the thermodynamic one. If the system 
crosses over smoothly from free subunits to disordered 
clusters, then there is no temperature at which assembly 
is efficient on the time scale i y i e id- 



B. Local equilibration 

We now return to the link between the kinetic crossover 
and the FDR. The general idea is that dynamics that is 
locally time-reversal symmetric allows disordered states 
to anneal into ordered states. This idea is not new (for 
example, see Ref. [22j |. especially its Fig. 1). However, 
the FDR provides a quantitative measure of this effect. 

The crystalline state of the two-dimensional system of 
discs is close-packed, with each particle bonded to six 
neighbors. During assembly, the fraction of such par- 
ticles in a given cluster provides a measure of its crys- 
tallinity. As clusters form, there are many possible states 
with low crystallinity, and the system tends to visit these 
states quite frequently. The effectiveness of assembly de- 
pends on whether these states are able to anneal into 
crystalline clusters. This annealing becomes more diffi- 
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cult as the disordered clusters grow. For example, an- 
nealing the disordered clusters of Fig. [2] into crystallites 
requires highly co-operative processes with large activa- 
tion energies, while annealing small disordered clusters 
requires less co-operativity. 

Our results indicate that near optimal assembly, large 
disordered clusters are avoided because the system re- 
mains locally equilibrated at each stage of the assembly 
process (although the system is globally out of equilib- 
rium). At any stage of assembly, there will be a set of 
likely states. The condition of local equilibration is that 
the relative probabilities of these likely states reflect their 
relative Boltzmann weights. If this condition holds, the 
system avoids the disordered states that are precursors 
to the large disordered clusters of Fig. [2l For example, 
small disordered clusters have smaller Boltzmann weights 
than crystalline clusters of the same size, so local equi- 
libration suppresses the disordered states. On the other 
hand, if disordered states are likely at any stage of assem- 
bly, this indicates that they are not being annealed into 
crystallites, and are likely to evolve into larger disordered 
clusters. 

To link this argument with the FDR, we first demon- 
strate a link between local equilibration and an approxi- 
mate time-reversal symmetry. We consider two states C 
and C that are both likely at a given stage of assembly. 
The rate with which the system makes transitions from 
C to C is 

1 (C^C , ,t)=W(C , \C)p(C,t) (7) 

where p(C, t) is the probability that the system is in state 
C at time t, and W(C'|C) is the rate for transitions to 
state C given that the system is initially in state C. [The 
rate W(C|C) depends only on the dynamical rules of the 
model, while the rate j(C — » C, t) depends also on the 
state of the system at time t] . 

For models that obey detailed balance, we have 

W(C'\C)exp(fiEc>) = W(C\C')ex P {l3E c ), (8) 

Further, if the system is locally equilibrated then we have 

P (C, t) exp(/3£c) « P(C, t) exp(/3Ec>). (9) 

where C and C are likely states at this time. Thus, the 
rates for forward and reverse transitions between these 
states are equal: 

7(C-»C',t)» 7 (C'->C,*). (10) 

This relation is an approximate time-reversal symmetry 
of the locally equilibrated state, which holds on time 
scales for which the set of likely states is not changing 
significantly. 

The extent to which this approximate time-reversal 
symmetry holds is correlated with the degree of local 
equilibration, and hence with the extent to which the 
system is discriminating between high-energy disordered 
states and low-energy ordered ones. By avoiding the 



high-energy disordered states, the locally equilibrated 
system tends to assemble effectively. 

To link this local equilibration with the FDR, we show 
in the appendix that, for systems obeying detailed bal- 
ance, deviations from FDT arise from differences between 
the probabilities of trajectories and their time-reversed 
counterparts, during the time between perturbation and 
measurement. The key result is (IA5|) . We conclude that 
the FDR is a probe of the degree to which the system 
obeys time-reversal symmetry between times t w and t, 
and hence of the degree of local equilibration. 

Thus, our results for both capsid and disc systems 
(Fig. [6]) are consistent both with our hypothesis that 
the system falls out of local equilibrium at the kinetic 
crossover, and with our interpretation of the FDR as 
a measure of local equilibration. The parametric plots 
summarize the important features of the correlation and 
response functions, in a single system-independent plot, 
in which time and energy scales are rescaled away. The 
qualitative similarities in the behavior of the FDR are all 
the more remarkable given the different dimensionalities 
of the two models that we consider, and the very different 
structures of their assembled states. 



C. Generic and non-generic features of the FDR 

While the behavior of both capsid and disc systems 
are both consistent with our analysis above, there are 
important differences between the two panels of Fig. [5] 
In particular, at the peak of the assembly curve, the re- 
sponse in the disc system is larger than the corresponding 
response in the capsid system. 

The reason for this difference can be inferred from the 
states shown in Fig. [T] and [2J In the disc system, the 
crystallinity of the product is rather low at all temper- 
atures. Even small clusters typically explore many dis- 
ordered states before they form locally crystalline struc- 
tures. The system needs to be very close to local equi- 
librium in order to ensure that the ordered structures 
can be discriminated from the large number of disordered 
states. Thus, assembly is effective only when the FDR is 
close to unity. On the other hand, the directional interac- 
tions in the model capsid system impose quite stringent 
constraints on the local structure of the growing clus- 
ter. This reduces the possibility for stable disordered 
clusters, and discriminating between ordered and disor- 
dered states is easier. Thus, the system still assembles 
effectively even when deviations from local equilibrium 
are quite significant, and assembly is still effective even 
when deviations from FDT are quite large. 

Taking account of these differences, we emphasize the 
main feature of Fig. O the FDR is large above the kinetic 
crossover, and small below it. We expect this behavior 
to be preserved as long three conditions are met. Firstly, 
the observables used to construct the FDR should couple 
to the processes by which metastable disordered states 
are annealed into ordered ones. For example, if we had 
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measured the FDR in the capsid system using the cap- 
somer positions in place of their energies, then diffusive 
processes would dominate both correlation and response 
functions, and this response is not sensitive to the extent 
to which the bonds in the system are locally equilibrated. 

Secondly, we require that the assembled state of the 
system minimizes the free energy both globally and lo- 
cally. Many biological systems are believed to have 
funnel-shaped energy landscapes consistent with this 
constraint (23|. The models presented in this article also 
have this property. We believe that satisfying this con- 
straint contributes quite generally to good assembly, and 
it is therefore practical to bear it in mind when designing 
self-assembling systems. Of course, systems that violate 
this constraint do exist. For example, in three dimen- 
sions, minimization of the free energy of small clusters 
of spherical particles lead to icosahedral structures [24[ , 
while the crystalline phase has a close-packed structure. 
It is therefore possible for these particles to assemble 
into icosahedral structures while always remaining locally 
equilibrated. The FDR would be close to unity, but the 
system would never visit the 'correctly assembled' close- 
packed structure. 

The third condition that is required to ensure useful- 
ness of the parametric FDR plot concerns the time t used 
to construct it. The behavior of Fig. [6] depends weakly 
on the value of t, but changing its order of magnitude 
will lead to different behavior. In particular, at very low 
temperatures and for large t, the capsid system shows 
an FDR close to unity. This occurs because the sys- 
tem is locally equilibrated over a particular set of disor- 
dered states. However, in this case, the system would 
not be locally equilibrated while the disordered clusters 
were forming, so that FDR on that time scale would have 
been smaller than unity. In other words, detection of 
the relevant deviations from local equilibrium requires a 
measurement on the time scales during which those de- 
viations occur. 

These three conditions show that the application of 
the FDR to self-assembling systems requires some con- 
sideration of the relevant observables and time scales. 
However, for the systems studied in this article, meet- 
ing these conditions does not require careful tuning of 
model parameters, but only the kind of qualitative anal- 
ysis discussed in this section. This represents evidence 
in favor of the applicability of these methods to other 
self- assembling systems. 



V. OUTLOOK 

The arguments of Section IIVI seem general, and rela- 
tively independent of details of the system. Further tests 
of the links between efficient assembly, local equilibra- 
tion, and FDRs would be valuable, especially if FDRs 
could be measured experimentally. In principle, FDRs 
can be obtained whenever conjugate correlation and re- 
sponse functions can be measured. Measuring fluctua- 



tions and responses of local quantities, such as the energy 
of a single subunit, requires a high degree of experimen- 
tal control, but methods do exist in some systems. For 
example, Wang et al [25[ recently measured an FDR in 
a three dimensional glassy colloidal system. The diffu- 
sive correlation function is conjugate to the response of 
a single particle to a force in that case. Applying similar 
methods to ordering processes of spheres or discs would 
be analogous to our studies of the sticky disc system. 

Turning to biological systems, it would be possible 
to measure the degree of kinetic frustration in the fold- 
ing of biomolecules, either computationally in more de- 
tailed capsid models, or in systems such as the trpzip 
peptide [20|, or experimentally, in RNA folding, by a 
generalization of the experiment of [27]. In this latter 
case, the conjugate variables of force and displacement 
are already measurable, although obtaining good statis- 
tics for the correlations and responses as a function of 
both t and t w might be challenging. Results obtained in 
this way would complement information about the non- 
equilibrium dynamics obtained from analysis of the work 
distribution [28|, [29(. For example, the thermodynamic 
definitions of reversible and irreversible work are linked 
to the idea that non-equilibrium processes can occur with 
or without local equilibration. By characterizing the ex- 
tent to which particular degrees of freedom are locally 
equilibrated on particular time scales, FDRs provide an- 
other link between these thermodynamic ideas and the 
statistical mechanics of non-equilibrium trajectories. 
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APPENDIX A: TIME REVERSIBILITY, AND 
THE FDR 

In this appendix, we briefly consider a general stochas- 
tic system evolving between times t\ and tf, and show how 
deviations from the predictions of FDT come from trajec- 
tories (histories) which occur with probabilities that are 
different from those of their time reversed counterparts. 

Consider a stochastic system evolving between times 
ti and tf. The energy of a configuration C during this 
time period is given by E(C) = E (C) — hA{C), where 
h is a field, A is an observable, and A(C) its value 
in configuration C. The stochastic dynamics obey de- 
tailed balance with respect to the Boltzmann distribu- 
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tion p e q(C) oc e P E \ C >. The response of observable B to 
the field h is 



X (t i ,t i ) = Y,B(C l ) 

C(t) 



dP[C{t);h] 



d{f3h) 



(Al) 



h=0 



where the sum is over trajectories (histories) of the sys- 
tem, which we indicate by the function C(t); the initial 
and final configurations of the trajectory are C\ and Cf 
respectively; and P[C(t); h] is the probability of a trajec- 
tory, which includes the probability of its initial condi- 
tion. 

The property of detailed balance dictates that 



log 



p[c(ty,h] 

Pi(Ci) 



log 



P[C(t);h] 
Pi(Ci) 



0h[A(Ct) - A^)} 

-0[E(Ct) - E(C$] 
(A2) 



where Pi(Cj) is the probability of the initial condition of 
the trajectory C(t), and C(t) is the time-reversed coun- 
terpart of C{t). That is, C(t) = TC(U + tf - t), where 
the operator T reverses all quantities that are odd under 
time reversal, such as momenta. To enforce time-reversal 
symmetry of the equilibrium state, we assume that the 
energy and its perturbation are time-reversal symmet- 
ric: E(C) = E(TC) and A(C) = A(TC). We also assume 
that B(C) = B(TC) for convenience, although the same 
analysis can also be carried out without this assumption, 
leading to an analogous result. 

Using (dP[C(t);h]/dh) 
P[C(t);h](d/dh)logP[C(t);h], we substitute flXg) 
into (|A1[) , and obtain 

X (tf,ii) = (B{tf)[A(ft) - A(jbi)]) + 



d 



J2B(C f )P[C(t);0] —- log P[C(t); h] 



C(t) 



h=0 

(A3) 



where we have used (•) = J2c(t) -PPWK')- The 
fluctuation-dissipation theorem states that the first two 
terms are equal at equilibrium, so we define A%(tf, U) — 
x(tf,ti) — (B(t{)[A(tf) — A(ti)]) in order to measure de- 
viations from FDT. 



To obtain an informative expression for Ax(tf,t{), we 
use conservation of probability to write 



J2B(Ci)PW);h]=J2 B ( C >^ 



(A4) 



c(t) 



Thus, the derivative of the left hand side of (| A4[) with 
respect to h is zero. Noting that B{d) — B(C{), we 
subtract this derivative from the right hand side of (j A3 jl , 
arriving at 



Ax(tf,ti) = ^B(Cf) 
c(t) 



d(fih) 



log P[C(t);h] 



h=0 



{P[C(t);0]-P[C(t);0]} (A5) 
The purpose of (|A5|) is to show that if all trajectories C(t) 
have the same probabilities as their time-reversed coun- 
terparts C(t), then the second term vanishes, and FDT 
applies. This condition holds exactly only at equilibrium, 
but if the dynamics of the system are close to local equi- 
librium between times t; and tf , then the relative weights 
of forward and reverse trajectories will be similar, and 
deviations from FDT will be small. 
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